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^ ' Abstract 



The reparametrization-invariant generating functional for the unitary and causal perturbation 
theory in general relativity in a finite space-time is obtained. The region of validity of the Faddeev- 
Popov-DeWitt functional is studied. It is shown that the invariant content of general relativity as 
a constrained system can be covered by two "equivalent" unconstrained systems: the "dynamic" 
(with "dynamic" evolution parameter as the metric scale factor) and "geometric" (given by the 
Levi-Civita type canonical transformation to the action-angle variables where the energy constraint 
converts into a new momentum). 

"Big Bang", the Hubble evolution, and creation of matter fields by the "geometric" vacuum 
are described by the inverted Levi-Civita transformation of the geomeric system into the dynamic 
one. The particular case of the Levi-Civita transformations are the Bogoliubov ones of the par- 
ticle variables (diagonalizing the dynamic Hamiltonian) to the quasiparticles (diagonalizing the 
equations of motion). The choice of initial conditions for the "Big Bang" in the form of the Bogoli- 
ubov (squeezed) vacuum reproduces all stages of the evolution of the Friedmann-Robertson- Walker 
universe in their conformal (Hoyle-Narlikar) versions. 
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*u ■ 1. Introduction 



"Big Bang" as the beginning of the Hubble evolution of a universe is described as a pure classical 
phenomenon on the basic of particular solutions of the Einstein equations in general relativity in the 
homogeneous approximation. A strange situation consists in that the highest level of the theory, i.e., 
the Faddeev-Popov-DeWitt generating functional for unitary S-matrix |3J |2), neglects the questions 
about the evolution of a universe which are in the competence of the simplest classical approximation. 
There is an opinion that the solution of the " Big Bang" problem in quantum theory goes beyond the 
scope of the unitary perturbation theory and even of general relativity. To answer these questions, we 
need a more general theory of the type of superstring ||. 

According to another point of view, the reason of theoretical difficulties in understanding the 
"Big Bang" phenomenon is not the Einstein theory, but the non-invariant method of its quantization. 
In particular, for the Faddeev-Popov-DeWitt unitary S-matrix the non-invariant coordinate time is 
considered as the time of evolution, whereas an observer in a universe can observe and measure only 



f 



invariants of group of diffeomorphisms of the Hamiltonian dynamics, which includes reparametrizations 
of the coordinate time [J|, [|, ^, . 

In the present paper we try to construct the unitary S-matrix for general relativity in a finite space- 
time in terms of the reparametrization-invariant evolution parameter, and to answer the questions: 
What do Quantum Universe and Quantum Gravity mean? What is the status of the "Big Bang" evolu- 
tion in quantum theory? What does creation of Quantum Universe mean? on the level of perturbation 
theory, using the scheme of the time-reparametrization- invariant Hamiltonian reduction M. 



In the context of the Dirac generalized Hamiltonian theory for constrained systems || 10, pf , 
this scheme means the explicit resolving of the first class constraints to determine the constraint-shell 
action directly in terms of invariants. In other words, we use the invariant reduction of the action (to 
obtain an equivalent unconstrained system) instead of the generally accepted non-invariant reduction 
of the phase space by fixing gauges 0, ||, [ll| (see Fig. 1). 

The example of the application of such an invariant reduction of the action is the Dirac formulation 
of QED |L3| directly in terms of the gauge-invariant (dressed) fields as the proof of the adequateness of 



the Coulomb gauge with the invariant content of classical equations. Recall that the invariant reduction 



of the action is the way to obtain the unconstrained Feynman integral [14] for the foundation of the 
intuitive Faddeev-Popov functional integral in gauges theories [O] and to reveal collective excitations 
of the gauge fields in the form of zero-modes of the first class constraints (see Fig. 2). 

A constructive idea of the considered invariant Hamiltonian reduction of general relativity is to 
introduce the dynamic evolution parameter as the zero-mode collective excitation of metric @, || |(| 0, 



16, 17, [18|, |19(| . This dynamic evolution parameter can be identified with the zero- Fourier harmonic 
of the space- metric determinant H || (treated in cosmology as the cosmic scale factor), whereas its 
conjugate momentum, i.e, the second (external) form, plays the role of the localizable Hamiltonian of 
evolution. 

The separation of this zero-mode evolution parameter on the level of the action allows us to 
determine also the invariant geometric time formed by averaging the time-like component of a metric 
over the space coordinates ]|, |(| . 

An observer reveals the evolution of the universe (with "Big Bang") as the dependence of the 
dynamic evolution parameter (i.e. cosmic scale factor) on the geometric time. 

The evolution of both a classical and a quantum universes in terms of the geometric time has 



the form of the canonical transformations [|20|, 21, 22 1 of the initial dynamic variables to a new set of 



variables for which the total energy constraint becomes a new momentum, and its conjugate variable 
(i.e., a new dynamic evolution parameter) coincides with the geometric time onto equations of motion. 

The content of the paper is the following. Section 2 is devoted to the reparametrization-invariant 
Hamiltonian reduction of general relativity. In section 3 we construct the generating functional for 
the unitary perturbation theory which includes "Big Bang" and Hubble evolution. In Section 4, "Big 
Bang" and Hubble evolution are reproduced in lowest order of perturbation theory. In Section 5, we 
research the conditions of validity of the conventional quantum field theory in the infinite space-time 
limit. 

2. Reparametrization-invariant Hamiltonian dynamics of GR 
2.1. Action and variables 

General relativity (GR) is given by the singular Einstein-Hilbert action with the matter fields 

V 2 d^„\ , r .1 ( ..2 _ 3 



W{g\n) = J d A x^-g[-^R(g) + C f ] (^fi 2 = M 



Planck 
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and by a measurable interval 

{dsf e = g aP dx a dx p . (2) 
They are invariant with respect to general coordinate transformations 



x 



x'n=x'(xo,x 1 ,x 2 ,x 3 ). (3) 



The generalized Hamiltonian approach to GR was formulated by Dirac and Arnovit, Deser and Misner 



[23 1 as a theory of system with constraints in 3 + 1 foliated space-time 

(dsf = g^dx^dx" = N 2 dt 2 -^g^dxi (dx* = dx { + N l dt) (4) 

with the lapse function N(t,x), three shift vectors N z (t,x), and six space components ^gij(t,x) 
depending on the coordinate time t and the space coordinates x. The Dirac-ADM parametrization 
of metric (||) characterizes a family of hypersurfaces t = const, with the unit normal vector u a = 
(1/iV, —N k /N) to a hypersurface and with the second (external) form 

-A^-A^) (5) 

that shows how this hypersurface is embedded into the four-dimensional space-time. 
Coordinate transformations conserving the family of hypersurfaces t = const. 

t^t = t(t); Xi -> Xi = Xi(t,x 1 , x 2 ,x 3 ) , (6) 

dt ~ u ^ ri dx k dt dx k dx l 

N = iV^; N k = N l - = - -=r (7) 

dt oxi dt oxi dt 

are called a kinemetric subgroup of the group of general coordinate transformations (|3|) [||, 24 ] . The 
group of kinemetric transformations is the group of diffeomorphisms of the generalized Hamiltonian 
dynamics. It includes reparametrizations of the nonobservable time coordinate t(t)) @ that play the 
principal role in the procedure of the reparametrization-invariant reduction discussed in the previous 
Sections. The main assertion of the invariant reduction is the following: the dynamic evolution 
parameter is not the coordinate but the variable with a negative contribution to the energy constraint. 
(Recall that this reduction is based on the explicit resolving of the global energy constraint with respect 
to the conjugate momentum of the dynamic evolution parameter to convert this momentum into the 
Hamiltonian of evolution of the reduced system.) 

A negative contribution to the energy constraint is given by the space-metric-determinant loga- 
rithm. Therefore, following papers [17, ||, 18, ||, 16, 19] we introduce an invariant evolution parameter 



ifo(t) as the zero Fourier harmonic component of this logarithm (treated, in cosmology, as the cosmic 
scale factor). This variable is distinguished in general relativity by the Lichnerowicz conformal-type 
transformation of field variables / with the conformal weight (n) 

(n )f - ={n)f (Mt)_y n t (8) 



where n = 2, 0, —3/2, —1 for the tensor, vector, spinor, and scalar fields, respectively, / is so-called 



conformal- invariant variable used in GR for the analysis of initial data [16, 25]. In particular, for 
metric we get 

9,u(t,x) = {^^-g^x) {dsf = (^p) 2 [NW - Wg^dxi] . (9) 
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As the zero Fourier harmonic is extracted from the space metric determinant logarithm, the space 
metric g~ij(t,x) should be defined in a class of nonzero harmonics 



d 3 x log \\gij(t,x)\ \ = 0. 



(10) 



The transformational properties of the curvature R(g) with respect to the transformations @ lead to 
the action (IIT) in the form |5| 



fa 



W(g\n) = W(g\ m ) -JdtJ d A x^-{^-). 
h V 



(11) 



This form define the global lapse function Nq as the average of the lapse function N in the metric g 
over the kinemetric invariant space volume 



N (t) 



V 



Jd*: 



N(t,x) 



9 = det(®g) , 



V = J d 3 x 



(12) 



V 



where Vq is a free parameter which in the perturbation theory has the meaning of a finite volume 
of the free coordinate space. The lapse function N(t, x) can be factorized into the global component 
No(t) and the local one M(t,x) 



N(t,x)g- 1/2 ■= N (t)M{t,x) := N q , 
where N fulfills normalization condition: 



(13) 



(14) 



that is imposed after the procedure of variation of action, to reproduce equations of motion of the 
initial theory. In the Dirac harmonical variables [p| chosen as 



q lk = gg lk , 



the metric takes the form 



(da) 



Mt? i 



12 ( N 2 dt 2 - qijdVdx^ , (q = det(q ij )). 



(15) 



(16) 



The Dirac-Bergmann version of action (11) in terms of the introduced above variables reads |], ||] 



fa 

W = J dt^L + -d t (P oVo ) 



(17) 



L 



J d3x [J2 P F F '- N l Vi 

V ^ F 



where 



P o 0o - 



p2 

4V K ^ ' 



F f 



(18) 
(19) 
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H(ip ) = J d s xMH(<p Q ) 
is the total Hamiltonian of the local degrees of freedom, 



(20) 



Ufa) = ^q lJ q kl [wji - mm] + {3) R(g) +H f , (21) 

and 

Pi = 2{V k (q kl Tru) - Vi(q kl 7r kl )) + P if (22) 

are the densities of energy and momentum and TCf,Pf are contributions of the matter fields. In the 
following, we call the set of the field variables F ( |l9|) with the dynamic evolution parameter tpo the 
field world space. The local part of the momentum of the space metric determinant 

7r(t,a:) :=g y 7ry (23) 

is given in the class of functions with the non-zero Fourier harmonics, so that 

J d 3 xn(t,x) = . (24) 

The geometric foundation of introducting the global variable @ in GR was given in |l9| as 
the assertion about the nonzero value of the second form in the whole space. This assertion (which 
contradicts the Dirac gauge tt = 0) follows from the global energy constraint, as, in the lowest order 
of the Dirac perturbation theory, positive contributions of particle-like excitations to the zero Fourier 
harmonic of the energy constraint can be compensated only by the nonzero value of the second form. 

The aim of this paper is to obtain the dynamic "equivalent" unconstrained system in the field 
world space (F, (fo) by explicit resolving the global energy constraint and to study the possibility of the 
Hamiltonian description of the "equivalent" unconstrained system in terms of the reparametrization- 
invariant evolution parameter T defined by equation NqcH = dT. 



2.2. Local constraints and equations of motion 

Following Dirac || we formulate generalized Hamiltonian dynamics for the considered system (0). It 
means the inclusion of momenta for M and iVj and appropriate terms with Lagrange multipliers 



t2 



W 



D 



dt\L u + -d t (P <po) \, L D = L+ d 6 x{P N N + P N *N l - \ Pm - X l P N i). (25) 



We can define extended Dirac Hamiltonian as 



+ J d 3 x(X°PM + X i P Ni ). 



(26) 



The equations obtained from variation of W D with respect to Lagrange multipliers are called first 
class primary constraints 

PM = 0, P N , = 0. (27) 

The condition of conservation of these constraints in time leads to the first class secondary constraints 

/ d^xNH 



{h d ,p m } = h 



VqAA 2 



0, 



{H D ,P N ,}=P i = 0. 



(28) 
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For completeness of the system we have to include set of secondary constraints. According Dirac we 
choose them in the form 

Af(t,x) = l; N i (t,x)=0; (29) 

7r(i,x) = 0; X j ■= di(q- 1/3 q ij ) = 0. (30) 

The equations of motion obtained for our system are 

dF _ dH(<p Q ) dP F _ dH(ipo) 

dT dP F ' AT dF ' 1 ) 

where H((po) is given by the equation (^), and we introduced reparametrization-invariant geometric 
time T 

N dt = dT . (32) 
2.3. Global constraints and equations of motion. 

The physical meaning of the geometric time T, the dynamic variable tpo and its momentum is given 
by the explicit resolving of the zero- Fourier harmonic of the energy constraint 

-JL + i% ) = 0. (33) 



SN (t) W 

This constraint has two solutions for the global momentum Pq: 



(P )± = ±2y/V H{<p ) = H* ± . (34) 
The equation of motion for this global momentum Pq in gauge (^9|) takes the form 

SW E (dtp\ (P )± , [-— ( , Jd 3 xH H(<p ) 



5P \dTj± 2V v^ u " ryru, Vq y Q 

The integral form of the last equation is 

T(<p l}V2 )= [ d^p- 1 ' 2 ^). (36) 



V>1 



Equation obtained by varying the action with respect to <po follows independently from the set of 
other constraints and equations of motion. 



Equations (35), (pq) in the homogeneous approximation of GR are the basis of observational 
cosmology where the geometric time is the conformal time connected with the world time Tf of the 
Friedmann cosmology by the relation 

= MTl dT (37) 

fi 

and the dependence of scale factor (dynamic evolution parameter ipo) on the geometric time T is 
treated as the evolution of the universe. In particular, equation ( p5| ) gives the relation between the 
present-day value of the dynamic evolution parameter <Pq(Tq) and cosmological observations, i.e., the 
density of matter p and the Hubble parameter 



n e hub = ^ = ^f => MTo)=(^) :=^o /4 ( 0.6 < (Ql /4 ) exp < 1.2 ). (38) 

</>0 ¥>0 \rihubJ 
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The dynamic evolution parameter as the cosmic scale factor and a value of its conjugate momentum 

are 



(i.e., a value of the dynamic Hamiltonian) as the density of matter (see equations (35), ( |36| ) 
objects of measurement in observational astrophysics and cosmology and numerous discussions about 
the Hubble parameter, dark matter, and hidden mass. 

Our aim is to find the equivalent unconstrained Hamiltonian system that describes evolution of 
the field world space (F, <po) in terms of the geometric time T. 



2.4. Equivalent Unconstrained Systems 

Assume that we can solve the constraint equations and pass to the reduced space of independent 
variables (F*,P F ). The explicit solution of the local and global constraints has two analytic branches 



with positive and negative values for scale factor momentum Pq (34). Therefore, inserting solutions of 
all constraints into the action we get two branches of the equivalent Dynamic Unconstrained System 
(DUS) 



W° us 



[F\<Po\ 



¥2 , r 



OF* 
dcpo 



(39) 



where ipo plays the role of evolution parameter and H± (defined by equation ( |34| ) plays the role of 
the evolution Hamiltonian in the reduced phase space of independent physical variables (F*,P F ) with 
equations of motion 

dF* dH* ± dP* F dH* ± 



dip dPy dipo dF* ' ^ 

The evolution of the field world space variables (F*,ipo) with respect to the geometric time T is not 
contained in DUS (|39|). This geometric time evolution is described by supplementary equation for 
nonphysical momentum Pq (|35| ) that follows from the initial extended system. 

To get an equivalent unconstrained system in terms of the geometric time (we call it the Geometric 



Unconstrained System (GUS)), we need the Levi-Civita canonical transformation (LC) [^. 21, [2^] of 
the field world phase space 

(F*,P F \ip ,P ) => (Fq,Pq\Q ,U ) (41) 

which converts the energy constraint (^3|) into the new momentum Hq (see the similar transformations 
for a relativistic particle in Appendix A). 

In terms of geometrical variables the action takes the form 



W G = J dt\ J d 3 xJ2 P G F G 



n Q + N Ii Q + ^-S LC 
dt 



(42) 



where S LC is generating function of LC transformations. Then the energy constraint and the supple- 
mentary equation for the new momentum take trivial form 



n n 



o 



sw 



Equations of motion are also trivial 







dT 



dQo 
dt 



dQo = dT. 



0, 



dF, 



G 



dT 



and their solutions are given by the initial data 

p* _ p* o 

G G ' 



F, 



G 



TP* 
^G ■ 



(43) 



(44) 



(45) 
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Substituting solutions of and ( f44| ) into the inverted Levi-Civita transformations 

F* = F*(Q ,U \Fq,Pq), <po = MQoM*G,Pg) (46) 

and similar for momenta, we get formal solutions of ( ^0|) and (|36| ) 

_F* = F*(T, 0\Fq°, Pq°), P f = P f (T,0\F^,P^°), <p = m (T, 0|i£°, Pq)- (47) 

We see that the geometric time evolution of the dynamic variables is absent in DUS. The geometric 
time evolution of the dynamic variables can be described in the form of the LC (inverted) canonical 
transformation of GUS into DUS ©, ©. 

To obtain the geometric time Hamiltonian evolution, it is sufficient to use the weak form of Levi- 
Civita-type transformations to GUS (F*,P F ) =>■ (F,P) with a new constraint 

Uo-H(Qo,F,P) =0. (48) 

We get the constraint shell action 

-H(T,F,P)\, (49) 

that allows us to choose the initial cosmological data with respect to the geometric time. 

The considered invariant reduction reveals the difference of reparametrization-invariant theory 
from the gauge-invariant theory: in gauge-invariant theory the superfluous (longitudinal) variables 
are completely excluded from the reduced system; whereas, in reparametrization-invariant theory the 
superfluous (longitudinal) variables leave the sector of the Dirac observables (i.e., the phase space 
(F*, Pp)) but not the sector of measurable quantities: superfluous (longitudinal) variables become the 
dynamic evolution parameter and dynamic Hamiltonian of the reduced theory. 



W 



GUS 



J cLt\ J d 3 x^P 



~dT 



2.5. Quantization and the arrow of the time 

In quantum theory of GR (like in quantum theories of a particle considered in Appendix A), we get 
two Schrodinger equations 

iJ_y±(F\<po, <pi) = FiM^O^o, n) (50) 
d(p 

with positive and negative eigenvalues of Po an d normalizable wave functions with the spectral series 
over quantum numbers Q 

y + (F\ip , n ) = J2 a q< F \Q >< Q\n, v>\ > e&o - vi) (5i) 

Q 

^-(F|v?o, Vl ) =J2 A Q< F \Q >*< Q\V0> Vl >* 0(^1 - V>o) , (52) 



where < F\Q > is the eigenfunction of the reduced energy (£ 

HU<Po) < P\Q >= ±#(Q. n) < F \Q > (53) 

<Q\<Po,<Pl >=exp[-i J dcpE(Q,<p)], <Q\<Po,<fi >*= exp[i J dipE(Q,<p)] . (54) 
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The coefficient Aq, in "secondary" quantization, can be treated as the operator of creation of a 
universe with positive energy; and the coefficient Aq, as the operator of annihilation of a universe also 
with positive energy. The "secondary" quantization means ^4q,^4q,] = 5q : ql The physical states 
of a quantum universe are formed by the action of these operators on the vacuum < 0|, |0 > in the 
form of out-state (\Q >= Aq\0 >) with positive "frequencies" and in-state (< Q\ =< 0|Aq) with 
negative "frequencies". This treatment means that positive frequencies propagate forward ((po > c/?i); 
and negative frequencies, backward (ip\ > ipo), so that the negative values of energy are excluded 
from the spectrum to provide the stability of the quantum system in quantum theory of GR. In other 
words, instead of changing the sign of energy, we change that of the dynamic evolution parameter, 
which leads to the causal Green function 

G c (Fi, yi\F 2 , cp 2 ) = G + (Fx,<pi\F 2 ,<p2)0(<P2 ~ <Pi) + G-(F U (pi\F 2 , <p 2 )0(<pi - y> 2 ) (55) 
where G + (Fi, <p\\F 2 , (p 2 ) = G-(F 2 ,ip 2 \Fi,ipi) is the "commutative" Green function 

G+{F 2 ,ip2\F\,<Pi) =< Ol^-^l^yi^+^lkl,^)^ > (56) 

For this causal convention, the geometric time (|36| ) is always positive in accordance with the equations 
of motion (35) 

•fiO 

(d^) =± ^P ^ r ±(Pi»Vo) = ± J dipp-^tcp) >0 . (57) 

Thus, the causal structure of the field world space immediatly leads to the arrow of the geometric 
time ( |57|) and the beginning of evolution of a universe with respect to the geometric time T = 0. 

As it was shown in [0], the way to obtain conserved integrals of motion in classical theory and 
quantum numbers Q in quantum theory is the Levi-Civita-type canonical transformation of the field 
world space (F,(po) to a geometric set of variables (V,Qo) with the condition that the geometric 
evolution parameter Qq coincides with the geometric time dT = dQo (see Fig. 3). 

3. Reparametrization invariant path integral 

Following Faddeev-Popov procedure we can write down the path integral for local fields of our theory 



using constraints and gauge conditions (|27| - |30|) : 

F 2 



^1(^1,^0,^0,^0) = J D(F,P f )A s A t exp{iW} , (58) 



Fx 

where 



t,x \i<k f 



are functional differentials for the metric fields (jr, q) and the matter fields (pf,f), 

A s =Y[5(V l ))5( X j )det{V l ,x ] }, (60) 



A t = H5(H((x))5(7r)det{H(<po) - P,tt}, L = IA^E^A\ ( 6 i) 
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are the F-P determinants, and 

t2 



w 



J dt\J <Px\Y,P f f\ - P tp - N 



ti ivb 



(62) 



is extended action of considered theory. 

By analogy with SR considered in Appendix A we define a commutative Green function as an 
integral over global fields (Po,ipo) and the average over reparametrization group parameter iV"o 

G+iFuwlFw) = J J] I diP0d ^ dN ° J Z local (n,^|Po^o,iVo), (63) 

where 

N = N/2tt6(0), S(0)=[dN . (64) 



The causal Green function in the world field space (-F, y?o) is defined as the sum 

G c (Fi,(pi\F2,tp 2 ) = G + (F 1 ,ip 1 \F 2 ,ip 2 )6(p 1 - <P2) + G + (F 2 ,ip 1 \F 2 ,pi)9((p2 - tpi). (65) 
This function will be considered as generating functional for the unitary S-matrix elements J26| 



'•P2 

S[ipi,ip 2 ] =< out ((p 2 )\T v exp {-i I d<p(H}) } \(ipi) in >, (66) 



'■p\ 



where T v is a symbol of ordering with respect to parameter (po, and < out (<p 2 )\, \(<f) in > are states of 
quantum Univers in the lowest order of the Dirac perturbation theory (M = 1; N k = 0; g y = 5ij+hfj), 
Hj is the interaction Hamiltonian 

Hf = H*-H$, H* = 2y/v H(<p), H* = 2^V H (<p) , (67) 

Hq is a sum of the Hamiltonians of "free" fields (gravitons, photons, massive vectors, and spinors) 
where all masses (including the Planck mass) are replaced by the dynamic evolution parameter ipQ @. 
For example for gravitons the "free" hamiltonian takes the form: 



H ( V0 ) = J d\ + f^h T f ) ; {hi = 0; 8,1% = 0). (68) 



In order to reproduce Faddeev-Popov integral for general relativity in infinite space-time Q , one should 
fix the dynamic evolution parameter at its present-day value ipo = \i (p8|), remove all the zero- mode 
dynamics Pq = (£q = 0, Nq = 1, and neglect the surface Newton term in the Hamiltonian. We get 



Z FP (F l7 F 2 ) = Z loca i(Fi, F 2 \P = 0, Wexp = »,N = 1), (69) 



or 



F 2 

Z FP (F 1 ,F 2 ) = J D(F,P f )A s A t exp{iW fp }, (70) 
Fx 



where 



+ °° / \ 2 

W fp = J dtj d 3 x r£ p FF- Hf p (iJ,) 1 , HfM = H(p) - f^-didj 

— nun \ F / 



q lJ , (71) 
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and 

At = I] S(H(p))5(jr)det{H(ji), tt}. (72) 



t,x 



The F-P integral ( |70| ) is considered as the generating functional for unitary perturbation theory in 
terms of S-matrix elements 

f 7 \ 

S[—oo\ + oo] =< out|Texp < — i / dtHj(u) > |in > . (73) 

Strictly speaking, the approximation (^) is not a correct procedure, as it breaks the reparametriza- 
tion-invariance. The region of validity of FP integral (|70| ) is discussed in next sections. 



4. Evolution of "Free" Quantum Universe 
4.1. Dynamic unconstrained system 



Possible states of a free quantum universe in 5-matrix (36) are determined by the lowest order of the 
Dirac perturbation theory given by the well-known system of "free" conformal fields (JsJ) , ( |68|) in a 
finite space-time volume [27, 0] 



*2 ,r 

W Q E = J dt( J d 3 xY,PFF 



PoW - N [- 



P, 



^ + H & )} + ±d (P <pa) 



(74) 



where Hq is a sum of the Hamiltonians of "free" fields (gravitons (|6q), photons, massive vectors, 
and spinors) where all masses (including the Planck mass) are replaced by the dynamic evolution 
parameter ipo M. The classical equations for the action d73) 



dF 
dT 



dHo 
dP F 



dP F dH 



dT dF 



P = ±2 y/VoIh := ^ 



(75) 



contain two invariant times: the geometric T and the dynamic (p^ connected by the geometro-dynamic 
(back-reaction) equation 



dT 



PO= Vo 



Solving the energy constraint we get the action for dynamic system 



W E (constraint) = W ( 



D 



J d^{ J d^xJ^PFd^F 



(76) 



(77) 



that has two branches for a universe with a positive energy (Pq > 0), and a universe with a negative 
energy (Pq < 0) . We interpret the branch with negative energy as an " anti- universe" which propagates 
backward ((p < 0) with positive energy to provide the stability of a quantum system. 

The content of matter in the universe is described by the number of particles Np k and their energy 
u>f{<Po, k) (which depends on the dynamic evolution parameter ipo and quantum numbers k, momenta, 
spins, etc.). Detected particles are defined as the field variables F = f 



m = E 



Cf((p )exp(ikiXi 



(a+(-k)+a f (k) 



(78) 



11 



which diagonalize the operator of the density of matter 



Po = £ ^y^N fik , N f (a) = \{a+a- f + aja+) 



f,k 



(79) 



We restrict ourselves to gravitons (f=h) Ch(^>o) = ^oV^i Uh{tpo,k) = \[W and massive vector 



particles (f=v) C v (ipo) = 1, ui v (<po,k) = \Jk 2 + y 2 ifQ, where y is the mass in terms of the Planck 
constant. 



4.2. Geometric unconstrained system 

The equations of motion (^|) in terms of a + , a~ @ are not diagonal 



dT 



X ■= iXa f = -H af Xa f , Xaj = \* ) ] H af 



a f 



-iAf , -u a . 



where nondiagonal terms Af=h v are proportional to the Hubble parameter ( |38l) 



A - ^° 
A f=h - , 

¥>0 



CO,. 



(80) 



(81) 



The "geometric system" (b + ,b) is determined by the transformation to the set of variables which 
diagonalize equations of motion (75) and determine a set of integrals of motion of equations ([75]) (as 
conserved numbers {Q}). 

To obtain integrals of motion and to choose initial conditions for the universe evolution we use the 
Bogoliubov transformations [^] and define " quasi-particles" 

b + = cosh(r)e~* 9 a + - i sinh(r)e* 9 a, b = cosh(r)e* 9 a + % sinh(r)e _;i9 a + , (82) 



or 



Xb 



OXa , 



which diagonalize the classical equations expressed in terms of " particles" (a + , a) , so that the " number 
of quasiparticles" is conserved 



d(b+b) 
dt 



b = exp(— i J dTuJb(T))bo. 
o 



(83) 



Functions r and 9 in (p2|), and the quasiparticle energy ujb in (^3() are determined by the equation of 
diagonalization 



i-^x b = Ho-^o - o- l H a d]xb = 



Xb 



V 0, -u b J 



in the form obtained in |?j 

Ofb = (ujj — 6'f) cosh(2r^) — (Af cos 26 f) sinh(2ry) 



(84) 



(85) 



= (ojf- 6»})sinh(2r / ) - (A/ cos 20/) cosh(2r/) 



-A / sin26' / . 
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Equations ([81])- (pq ) are closed by the definition of " observable particles" in terms of quasiparticles 

X>/0){a| «/} 
H f r + i ri+i. i t n * 



p(<p) = _Ji = i ; { a + a} = {6,[6o} cosh 2r - -(6 +2 - b 2 ) sinh 2r (86) 

K VQ 2 



with 



s> /6 = y/{ Uf - e' f y + - a}, e' f = -l(x- f ) t 1 " ^Af rV "' cosh ( 2r /) = • ( 8? ) 

The constrained system in terms of geometric variables is described by the action 
W G 



L i_iJ^_ -1/2 cosh(2r f ) = ^ £ 



-Db + 5>£(Q )iV>(&) 
/ 



where the new dynamic evolution parameter Qo coincides with geometric time T on the equations of 
motion 

fg!=0 - dOo = <iT - (89) 

olio "I 

Reduction of this system leads to the weak version of Geometric Unconstrained System ( |49| ) 



W™* = JdT ~ (bd T b + - b + d T b) f -^J b (T)N f (b)\ . (90) 

We choose the initial data appropriate for the dynamics described by GUS 
4.3. Quantization 

The initial data boj&o of quasiparticle variables ( |B3| ) form the set of quantum numbers in quantum 
theory. 

Let us suppose that we manage to solve equations (|S3"1)- ( p7|) with respect to the geometric time T 
in terms of conserved numbers bt , bo. This means that the wave function of a quantum universe can 
be represented in the form of a series over the conserved quantum numbers Q = rif^ =< Q\btbf\Q > 
of the Bogoliubov states 



{/■ ] (b + ) n 

-i / dTn f uj b {T) \ K -^= 
o J VUf 



|0> 6 • (91) 



In this geometric system, we have an arrow of the geometric time T for the universe 

T + (<f>2,<pi) = f d<pp(^r 1/2 > 0, p 2 > (pi , (92) 



and for anti-universe 



91 <pi 

T-(p 2 ,Pi) = - I dpp(p)- 1 / 2 = f dpp^Y 1 ' 2 > 0, tp!> (p2. (93) 



<Pl ¥2 
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The dynamic system ( |77[ ) of particle variables a + ,a is connected with the geometric one by the 
Bogoliubov transformations. Using these transformations we can find wave functions of a universe, 
for if2 > (fi and an anti-Universe, for ip\ > if2 

* Q (T) = A+^T+fa, ipiWfa -ipi) + AQ$* Q (T-{ip 2 , ipi))e(<pi - if*) , (94) 

where the first term and the second one are positive (Po > 0) and negative (Pq < 0) frequency parts of 
the solutions with the spectrum of quasiparticles Aq is the operator of creation of a universe with a 
positive "frequency" (which propagates in the positive direction of the dynamic evolution parameter) 
and Aq is the operator of annihilation of a universe (or creation of an anti-universe) with a negative 
"frequency" (which propagates in the negative direction of the dynamic evolution parameter). 

We can see that the creation of the universe in the field world space and the creation of dynamic 
particles by the geometric vacuum (6 + |0 >= 0) are two different effects. 

The second effect disappears if we neglect gravitons and massive fields. In this case, dp/ dip = 0, 
and one can represent a wave function of Universe in the form of the spectral series over eigenvalues 
Pq of the density p 

*(/l¥>2,¥>i)= (95) 



E 



exp \-i(<p 2 - W ) £ < f\Q > + exp Ufa - ^ £ S^l < f\Q >* 



where < f\Q > is a product of normalizable Hermite polynomials. 



4.4. Evolution of quantum universe 

The equations of diagonalization (|84|) for the Bogoliubov coefficients (p2|) and the quasiparticle energy 
cjf, ( p5[) play the role of the equations of state of the field matter in the universe. We can show that 
the choice of initial conditions for the "Big Bang" in the form of the Bogoliubov (squeezed) vacuum 
6|0 >b= reproduces all stages of the evolution of the Friedmann-Robertson- Walker universe in their 
conformal versions: anisotropic, inflation, radiation, and dust. 

The squeezed vacuum (i.e., the vacuum of quasiparticles) is the state of "nothing". For small tp 
and a large Hubble parameter, at the beginning of the universe, the state of vacuum of quasiparticles 
leads to the density of matter @ 

b < P (a\ «)> b =^(^ + ^), B = \ (96) 

where </?(0) is the initial value, and po is the density of the Casimir energy of vacuum of " quasiparticles" . 
The first term corresponds to the conformal version of the rigid state equation (in accordance with the 
classification of the standard cosmology) which describes the Kasner anisotropic stage T±(ip) ~ ±(/? 2 
(considered on the quantum level by Misner |2^] ) . The second term of the squeezed vacuum density (p6|) 
(for an admissible positive branch) leads to the stage with inflation of the dynamic evolution parameter 
<p with respect to the geometric time T 

<p(T) (+) ~ p(0)exp[Tv^M0)]. 

It is the stage of intensive creation of "measurable particles" . After the inflation, the Hubble parameter 
goes to zero, and gravitons convert into photon-like oscillator excitations with the conserved number 
of particles. 
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At the present-day stage, the Bogoliubov quasiparticles coincide with particles, so that the mea- 
surable density of energy of matter in the universe is a sum of relativistic energies of all particles 

po^ = v Q =Y,y^ + y> 2 (n (97) 

rtf 

where yy is the mass of a particle in units of the Planck mass. The case of massless particles (y = 
0,po(tp) = constant) correspond to the conformal version of radiation stage of the standard FRW- 
cosmology. And the massive particles at rest (k = 0, po(^f) = PbarionsV 9 // 1 ) corresponds to the conformal 
version of the dust universe of the standard cosmology with the Hubble law 

tp=±^^ L P±\T) = y — j T ; q =— T = -. (98) 

The dynamic evolution parameter is expressed through the geometric time of a quantum asymptotic 
state of the universe \out > and conserved quantum numbers of this state: energy E ou t and density 
Po = Eout/Vo. 

It is well-known that E ov ± is a tremendous energy (10 79 GeU) in comparison with possible real and 
virtual deviations of the free Hamiltonian in the laboratory processes: 

H = E out + 5H , < out\6H \in >« E out . (99) 

We have seen that the dependence of the scale factor tpo on the geometric time T (or the "relation" of 
two classical unconstrained systems: dynamic and geometric) describes the "Big Bang" and evolution 
of a universe. 

Therefore, from the point of view of unconstrained system "Big Bang" is effect of evolution of the 
geometric interval with respect to the dynamic evolution parameter which goes beyond the scope of 
Hamiltonian description of a single classical unconstrained system. 

Reparametrization-invariant dynamic of General Relativity is covered by Geometric and Dynamic 
Unconstrained Systems connected by the Levi-Civita transformation of fields of matter into the vacuum 
fields of initial data with respect to geometric time (see Fig. 3). 



5. QFT limit of Quantum Gravity 

The simplest way to determine the QFT limit of Quantum Gravity and to find the region of validity of 
the FP-integral(^) is to use the quantum field version of the reparametrization-invariant integral (|63l) 
in the form of S-matrix elements [^] (see (|66|), (|67[)). We consider the infinite volume limit of the 
S-matrix element (p7|) in terms of the geometric time T for the present-day stage T = Tq, <p(Tq) = u, 
and T{ipi) = Tq — AT,T((p2) =Tq + AT = T OU f. One can express this matrix element in terms of the 
time measured by an observer of an out-state with a tremendous number of particles in the Universe 
using equation (||) dip = dT outy J p ou t and approximation ( p9| ) to neglect "back-reaction". In the 
infinite volume limit, we get from 



dpo[#j] = 2cVo (y^ (#o + Hj) - VVoH^j = dT out [FHj + 0{l/E out )} (100) 
where Hj is the interaction Hamiltonian in GR, and 

^i/W^idrk <101) 
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is a multiplier which plays the role of a form factor for physical processes observed in the "laboratory" 
conditions when the cosmic energy E ou t is much greater than the deviation of the free energy 

5H = H - E out ; (102) 

due to creation and annihilation of real and virtual particles in the laboratory experiments. 

The measurable time of the laboratory experiments T2 — T\ is much smaller than the age of the 
universe To, but it is much greater than the reverse "laboratory" energy 5, so that the limit 

T{(p2) +00 

dT OVLt => J dT out 



is valid. If we neglect the form factor (101) that removes a set of ultraviolet divergences, we get 
the matrix element ( |73|) that corresponds to the standard FP functional integral ( |70|) and S'-matrix 
element ( |73| ) where the coordinate time is replaced by the geometric (conformal) time t — ► T out : 

S[— 00 1 + 00] =< oui|Texp |— i J dT out FHi([i) | |in > ^ F = l) . (103) 

Thus, the standard FP-integral (i.e., the Hamiltonian description of the evolution of fields with respect 
to the geometric time) appears as the nonrelativistic approximation of tremendous mass of a universe 
and its very large life-time (see Fig. 4). Now, it is evident that the conventional FP unitary S-matrix 
are not valid for the description of the early universe given in the finite spatial volume and the finite 
positive interval of geometrical time when the early universe only begins to create matter. 

On the other hand, we revealed that standard quantum field theory (103) is expressed in terms 



of the conformal- invariant Lichnerowicz variables and coordinates including the conformal time (T ou t) 
as the time of evolution of these variables. The reparametrization-invariant description of the "Big 
Bang" evolution distinguishes conformal variables (||) and coordinates. The conformal invariance of 
the variables can testify to the conformal invariance of the initial theory of gravity. Such the theory 
can be a scalar version of the Weyl conformal invariant theory j|, [f| ^0| (that dynamically coincides 
with the Einstein General Relativity) where these conformal variables, coordinates, geometric time 
T, and the conformal Hubble parameter TL c hub = f' /f can be considered as measurable quantities. In 
particular, to get "accelerating" universe with q > in the dust stage (see (|98|)), is enough to count 
that we measure the relative Weyl time (T ou t) of quantum field theory ( |103| ). 



6. Conclusions 

The main result of the paper is the reparametrization-invariant generating functional for the unitary 
and causal perturbation theory in general relativity in a finite space-time. We show that the classical 
cosmology of a universe and the Faddeev-Popov-DeWitt functional correspond to different orders of 
decomposition of this functional over the inverse "mass" of a universe. 

This result is based on the assertion that the measurable time in any reparametrization-invariant 
system is not the coordinate time, but the time-like dynamic variable of an extended phase space of 
this system (of the type of the conformal scale factor). Accordingly, the measurable Hamiltonian is a 
solution of the energy constraint with respect to the conjugate momentum of this dynamic evolution 
parameter |]|, |6| . This definition of the dynamic evolution parameter and Hamiltonian supposes the 
reduction of an action for constructing an "equivalent" unconstrained system. 
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The second assertion is that such an unconstrained dynamic system cannot cover the physical 
content of a relativistic reparametrization-invariant system. This content can be covered by two 
"equivalent" unconstrained systems - the mentioned above dynamic system and the geometric sys- 
tem constructed by the Levi-Civita type canonical transformation so that a new dynamic evolution 
parameter coincides with the geometric time. 

The " geometric" variables are the Bogoliubov quasiparticles which diagonalize equations of motion 
and give cosmological initial conditions. The coefficients of the Bogoliubov transformations, in the 
conventional QFT perturbation theory, determine the evolution density of matter. For the vacuum 
initial data this evolution reproduces all stages of the standard FRW evolution of the universe in their 
conformal versions. 

Consistent QFT limits of the generating functional in classical gravitation and cosmology and 
the conventional quantum field theory (in the form of the Faddeev-Popov generating functional for 
a infinite space-time) can be fulfilled for a tremendous value of the universe mass and the universe 
life-time (see Fig. 4). The quantum field theory limit distinguish the conformal treatment of general 
relativity developed in || |(| [3(]] . 
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APPENDIX A. Special Relativity 

To answer the question: Why is the reparametrization-invariant Hamiltonian reduction needed?, 
let us consider relativistic mechanics || in the Hamiltonian form 

T2 

r N 

W[P,X\N\n,r 2 ] = J dT[-P fl X»-—(-PZ + m 2 )} . (Al) 

Tl 

This action is invariant with respect to reparametrizations of the coordinate evolution parameter 

t^t' = r(r), N'dr 1 = NdT (A.2) 
given in the one-dimensional space with the invariant interval 

T 

dT := NdT, T = J dfN(f) . (A.3) 

o 

We called this invariant interval the geometric time []|] whereas the dynamic variable Xq (with a 
negative contribution in the constraint) we called dynamic evolution parameter. 

In terms of the geometric time ( |A.3|) the classical equations of the generalized Hamiltonian sys- 
tem ( |A.1| ) takes the form 

^i = 5i, ^ = 0, P 2 -m 2 = 0. (A.4) 

dT m dT ' M y ' 
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The problem is to obtain the equivalent unconstrained theories directly in terms of the invariant 
times Xq or T with the invariant Hamiltonians of evolution with respect to these times. The solu- 
tion of this problem is called the dynamic (for Xq), or geometric (for T) reparametrization-invariant 
Hamiltonian reductions. 

The dynamic reduction of the extended system ( |A,1| ) means the substitution, into it, of the explicit 
resolving of the energy constraint (-P 2 + m 2 ) = with respect to the momentum Pq with a negative 
contribution 

SW 



.... => Po = ±V m2 + ^ 2 - ( A - 5 ) 



In accordance with the two signs of the solution ( A, 5 ), after the substitution of ( A. 5 ) into ( A.l ), we 



have two branches of the dynamic unconstrained system 

X (t 2 )=X (2) 

W (constraint) = W±[Pi,Xi\X (l), X (2)] = J dX 

X (tl)=Xo(1) 



P^~ T M'P + P? 

dXn 



(A.6) 



The role of the time of evolution, in this action, is played by the variable Xq which abandons the Dirac 
sector of " observables" Pi,Xi, but not the sector of "measurable" quantities. At the same time, its 
conjugate momentum Pq converts into the corresponding Hamiltonian of evolution, values of which 
are the energy of a particle. 

This invariant reduction of the action gives the "equivalent" unconstrained system together with 
definition of the invariant evolution parameter (Xq) corresponding to a non-zero Hamiltonian Pq. 

Thus, we need the reparametrization-invariant Hamiltonian reduction to determine the invariant 
evolution parameter and its invariant Hamiltonian for reparametrization-invariant systems. 

In quantum relativistic theory, we get two Schrodinger equations 

l J r * {±) (X\P) = ±Jm* + P?* i±) (X\P) , (A.7) 
with positive and negative values of Po an d normalized wave functions 



± 



(X\P) = £ P 5t%- exp(-iP^), ([A P , A+) = 5 3 (Pi ~ Pi)) ■ (A.« 



The coefficient Ap~, in the secondary quantization, is treated as the operator of creation of a particle 
with positive energy; and the coefficient Ap, as the operator of annihilation of a particle also with 
positive energy. The physical states are formed by action of these operators on the vacuum < 0|, |0 > 
in the form of out-state ( \P >= ^4p|0 > ) with positive frequencies and in-state ( < P\ =< 0\Ap ) with 
negative frequencies. This treatment means that positive frequencies propagate forward (X02 > -^01); 
and negative frequencies, backward (Xqi > XQ2), so that the negative values of energy are excluded 
from the spectrum to provide the stability of the quantum system in QFT [26]. For this causal 



convention the geometric time ( |A.3| ) is always positive in accordance with the equations of motion 

n r P \ 77} 777 

T(X 02 ,X 01 )=± , (X 02 -X 01 )>0 (A.c 



dX J± Jp2 + m 2 ' /p2 +m 2 



In other words, instead of changing the sign of energy, we change that of the dynamic evolution 
parameter, which leads to the arrow of the geometric time (A. 9) and to the causal Green function 



G c (X)=G + (X)6(X ) + G^(X)e(-X ) = iJ -^exp(-iPl) p2 _^_, e , (A.10) 
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where G+(X) = G~(—X) is the "commutative" Green function [26] 

G+(X) = I -f-^exp(-iPX)5(P 2 -m 2 )9(P ) = — [ d 3 Pd 3 P' < 0|*_(X|P)* + (0|P / )|0 > . 
J (2tt) 6 2it J 

(A.ll) 

To obtain the reparametrization-invariant form of the functional integral adequate to the considered 
gauge-less reduction ( |A.6| ) and the causal Green function ( [A. 10 ), we use the version of composition 
law for the commutative Green function with the integration over the whole measurable sector X\^ 



G + (X -X )= / G + (X - X^G+iX! - X )dX 1 



G A 



2vrJ(0) ; ' 



(A.12) 



where 5(0) = J dN is the infinite volume of the group of reparametrizations of the coordinate r. Using 
the composition law n-times, we got the multiple integral 



/n 
G + (X -X 1 )l[G + (X k -X k+1 )dX k , (X n+1 =X ). (A.13) 
fc=i 

The continual limit of the multiple integral with the integral representation for ^-function 

5(P 2 - m 2 ) = — f dNexp[iN(P 2 - m 2 )} 
2tt J 

can be defined as the path integral in the form of the average over the group of reparametrizations 

a+(x) = x y m *gpw> j |^ (T) n ( smi ) ] eMmP , x|wln , 



x(n)=o 



(A.14) 

where N = N/2tt8(0), and W is the initial extended action ( |A.1| ), 

This functional integral has the form of the average over the group of reparametrization of the 
integral over the sector of "measurable" variables P^^X^. 

The Hamiltonian unconstrained system in terms of the geometric time T can be obtained by the 
canonical Levi-Civita - type transformation p0| , gj , |22j 



(-f/Ltj X^) =>• (U^Q^j 



(A.15) 



to the variables (n^, Q^) for which one of equations identifies Qq with the geometric time T. This 
transformation [2C] converts the constraint into a new momentum 



1 



in 



Pi 



no — - — [Pq — P{ ], rij — Pi, Qo — Xq—, Qi — Xi — x — 

2m P P 



and has the inverted form 



±\2mU + W, P l =U l , X = ±Q C 



2mU + U 2 



m 



Xi — Qi + Qo - 



n, 



■/?•) 



After transformation ( |A.16| ) the action ( [A.lj) takes the form 

TO 

d 



W 



dr 



-IW 1 - N(-IL + -) - — s 



Ic 



dr 



s lc = (Q n ). 



(A.16) 



(A.17) 



(A.18) 



n 
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The invariant reduction is the resolving of the constraint Ho = m/2 which determines a new Hamilto- 
nian of evolution with respect to the new dynamic evolution parameter Qq, whereas the equation of 
motion for this momentum Uq identifies the dynamic evolution parameter Qq with the geometric time 



T (Qo = T). The substitution of these solutions into the action ( A.18| ) leads to the reduced action of 
a geometric unconstrained system 

W(constraint) = W G [U,Q i \T 1 ,T2] = J dT (ll^ - ^ - (S lc = Q„y), (A.19) 

Ti 

where variables IT , Qi are cyclic ones and have the meaning of initial conditions in the comoving frame 
SW dQi „ „( ) 5W dUi „ , 0) , A . 

The substitution of all geometric solutions 

q = t, n = ^, H i = lf ) =P i , Qi = Qi° ] (A.21) 



into the inverted Levi-Civita transformation ( |A.17 ) leads to the conventional relativistic solution for 
the dynamical system 

r(°) v (r\ - rr P ° v fr\ _ , -r Pi 



P = ±V m2 + i f' P i = U i X (T)=T^, X i (T)=X\" ) +TA (A.22) 
v m m 

The Schrodinger equation for the wave function 

A^c^q.ijj.) = ^ylc^QilILi), ^(T.Qillli) =exp(-iT^)exp(mf ) Q 4 ) (A.23) 

contains only one eigenvalue m/2 degenerated with respect to the cyclic momentum rij. We see 
that there are differences between the dynamic and geometric descriptions. The dynamic evolution 
parameter is given in the whole region — oo < Xq < +00, whereas the geometric one is only positive 
< T < +00, as it follows from the properties of the causal Green function ( A.10| ) after the Levi-Civita 
transformation ( A.16| ) 



-co 



-co 



Two solutions of the constraint (a particle and antiparticle) in the dynamic system correspond to a 
single solution in the geometric system. 

Thus, the reparametrization-invariant content of the equations of motion of a relativistic particle 
in terms of the geometric time is covered by two "equivalent" unconstrained systems: the dynamic and 
geometric. In both the systems, the invariant times are not the coordinate evolution parameter, but 
variables with the negative contribution into the energy constraint. The Hamiltonian description of a 
relativistic particle in terms of the geometric time can be achieved by the Levi-Civita-type canonical 
transformation, so that the energy constraint converts into a new momentum. Whereas, the dynamic 
unconstrained system is suit for the secondary quantization and the derivation of the causal Green 
function that determine the arrow of the geometric time. 
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FIGURE CAPTIONS 

Fig. 1. The tree ot modern theoretical physics grew from two different roots ("particle" and 
"field") which gave the VARIATIONAL method and SYMMETRY principles for formulating modern 
physical theories as constraned systems. To obtain unambigious physical results, one should construct 
Equivalent Unconstrained Systems compatible with the symplest variational method. It is just the 
problem discussed in the present paper. 

Fig. 2. An equivalent unconstrained system W*(p*,q*) can be obtained in the case when the 
operations of varying and constraining commute with each other. The next problem is to establish 
the range of validity of the standard Faddeev-Popov (FP) integral. 

Fig. 3. Reparametrization-invariant dynamics of General Relativity is covered by the Dynamic 
Unconstrained Systems (DUS) and the Geometric Unconstrained Systems (GUS) connected by the 
Levi-Civita transformation of fields of MATTER into the vacuum fields of initial data with respect to 
geometric TIME. 

Fig. 4. Reparametrization-invariant dynamics of General Relativity. The Big Bang of Quantum 
Universe from point of view of Geometric and Dynamic Unconstrained Systems connected by Levi- 
Civita canonical transformations. 



22 




Fig. 1 



Y-M 

CONSTRAINED 
SYSTEM 



CONSTRAINT-SHELL 



UNCONSTRAINED 
SYSTEM 




Z Fp = J dpdqA Fp 8 (gauge)c) (const)£ 



iW 



sing 



UNITARY PERTURBATION THEORY 



Fig. 2 




Fig. 3 




LABORATORY EXPERIMENT 



Fig. 4 



